ISSN 2319 - 6629 



Volume 5, No.3, April - May 2016 

International Journal of Wireless Communications and Networking Technologies 

Available Online at http://warse.org/IJWCNT/static/pdf/file/ijwcnt02532016.pdf 

Availability of Free-Space Optical 
Systems Depending on Atmospheric Conditions 
and System Parameters 


Tsvetan Mitsev 1 , Yordan Kovachev 2 

faculty of Telecommunications at Technical University of Sofia, Bulgaria, Email: mitzev@tu-sofia.bg 
2 Faculty of Telecommunications at Technical University of Sofia, Bulgaria, Email: dakatapz@gmail.com 


ABSTRACT 

In this paper the availability of Free-Space Optical (FSO) 
systems is studied. In order to create mathematical model 
and calculate systems availability, meteorological data from 
different geographic regions is used to obtain statistical 
model of the atmospheric visibility S M - Utilizing the 
cumulative distribution functions ( CDF) of S M and the 
statistical distribution of the misalignments (Ap) between 
the laser beam axis and the center of the receiving antenna, 
a method for calculating FSO availability is proposed. 
Based on the derived analytical expressions, numerical 
results are obtained and FSO availability is analyzed. 
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1. INTRODUCTION 

The commercialization of FSO systems in the past years 
necessitates the improvement of their channel capacity and 
reliable work (availability). This poses the question for 
selecting the optimal parameters of the free-space optical 
communication system [1]. One of the causes for the 
decrease of the FSO reliability is the random shifts of the 
laser beam axis from its original direction. These 
misalignments can be caused by mechanical vibrations of 
the transmitting aperture [2, 3], strong atmospheric 
turbulence [4] or pointing errors [8]. They can significantly 
reduce the quality of the transmitted data and the FSO 
availability. Another factor for unreliable work of the FSO 
system is the atmospheric visibility S M [9]. Significant 
decrease of S M can cause system outage. 

The remainder of this paper is paper is organized as 
follows: Section 2 describes the factors that may cause 
random fluctuations in the laser beam direction; Section 3 
presents the problem, which is solved in this study. In 
Section 4 the models for calculating the minimal optical 
power in the photodetector, the optimal beam divergence 
angle (Section 4. 1), the atmospheric visibility (Section 4.2), 
the random vibration in the optical beam direction (Section 
4.3) and the availability (Section 4.4) are presented. Section 
5 contains the numerical simulations and Section 6 is the 
conclusion. 


2. REASONS FOR RANDOM FLUCTUATIONS OF 
THE LASER BEAM DIRECTION 


Some of the main reasons for linear displacements (Ap) of 
the optical beam from its original direction are: 

- Mechanical vibrations of the transmitter, caused by 
instability of the foundation, where FSO transceiver is 
mounted [13, 14] (FSO is located on vibrating platforms 
within industrial warehouses, near production sites, etc.) 

- Building sway (when the FSO system is located on top of 
very high buildings) 

- Turbulent eddies [8, 12] 

The first two cases are identical. Mechanical vibrations in 
the transmitter’s side cause a linear shift in the direction of 
the optical beam. The jitter in the direction of the laser 
beam propagation triggered by atmospheric turbulence is 
more complex physical process [4]. It is caused by the 
largest turbulence eddies, which sizes / ~ L 0 (/ 0 is the inner 
scale size of the turbulent eddies and Lo is the outer scale 
size) are proportional to the distance between the 
transmitter and the receiver. 



Figure 1: Change of the laser beam direction caused by turbulent 
eddies 


However, no matter the cause of the misalignment, it results 
in a linear shift Ap, between the center of the receiving 
aperture and the optical beam axis (Figure 1). So it can be 
generalized as a jitter in the initial direction of the laser 
beam propagation. 


When the transmitting aperture (TA) is shifted in any 
direction with a random angle y (Figure 2), the receiving 
aperture (RA), with radius R r , is displaced from the axis z of 
the optical beam. This misalignment Ap is also random, 
because, when z = Z: 

Ap = yZ (1) 
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Figure 2: Misalignment between the center of the receiving 
aperture and optical beam axis 


In Figure 2 p z is the radius of the optical beam at the plane 
of the receiver, which is defined by 7(0, z)/e 2 . Optical 
radiation intensity along the beam axis is 7(0, z). 



Figure 3: Distribution of the optical intensity at the plane of the 
receiver 


3. DEFINING THE PROBLEM 

Commercially available systems work with collimated 
optical beams. As a result there is a very high reserve of 
optical energy at the receiving end of the system, which 
leads to FSO working with very low bit error rates (BER < 
10 ' 20 ). Afore mentioned conditions (see Section 2) can 
cause misalignments between the laser beam axis and the 
center of the receiving aperture. This means that working 
with collimated laser beam increases the probability for 
system outage, caused by unbearable jitter in the direction 
of the laser radiation. The great reserve of optical energy, 
when using non-divergent optical beams provides the 
possibility to enlarge the divergence angle of the optical 
beam and therefore its radius p z at the plane of the receiving 
aperture [7, 10]. The larger beam radius at the receiver’s 
side allows for greater linear misalignments between the 
center of the receiving aperture and the optical radiation 
axis. This lowers of fully eliminates the probability of 
outage caused by random fluctuations in the initial 
direction of the laser beam. The FSO parameters, which can 
be manipulated until achieving their optimal values are: 
beam divergence angle, and respectively, laser beam radius. 

Consider the case shown in Figure 3. It depicts the radial 
distribution of the optical radiation intensity in the plane of 
the receiver’s aperture 7(p, z). The intensity 7^ corresponds 
to the minimum optical power through the aperture of the 
receiver, for which FSO system will work reliably. The 
respective magnitude of the beam radius p max defines the 
angle 0 max . This is the value of bearable angular 
misalignment (the respective maximal linear misalignment 
is Pmax) of the laser beam axis from its main direction (0 = 0) 
due to different random factors. The divergence angle of the 
optical radiation is represented with 0 t . 


In Figure 3 two distributions of 7(p, z) are represented with 
7i(p, z) and 7 2 (p, z), respectively. They correspond to two 
divergence angles 0 t ,i and 0 t ,2 (0 t ,i < 0t,2) of the transmitted 
laser beam. The figure outlines the case of 0 m ax(0t,i) > 
0max(0t,2), but it is obvious that this trend will continue to the 
limit value 0 t;Opt . After this point, increasing 0 t will result in 
decrement in the magnitude of 0 max . The task of the 
following analysis is to calculate 0 t;Opt , and to investigate the 
availability of FSO system working with this optimal value 
of the beam divergence angle. 

4. MATHEMATICAL DESCRIPTION 


4.1. Calculation of the Minimal Optical Power in the 
Photodetector and the Optimal Beam Divergence Angle 


The filed distribution at receiver’s site, as shown in Figure 
3, is Gaussian [5, 11], therefore 


l(p,z)= /(0,z)exp 


( 






( 2 ) 


When 0 t is optimal, the effective radius is equal to the 
maximal linear misalignments (p = p max ) and field intensity 
is minimal 7 = 7 min , i.e. 

f 2 

r max 


= /(0,z)exp 


-2 


pIU) 


(3) 


Optical radiation intensity along the axis of the laser beam 
depends on the characteristics of the transmitter and the 
atmospheric communication channel [1]: 

2M ( 4 ) 

TT.PzUj 


/(0,z) = - 


where x t represents the losses in the receiving antenna and 
x a is the transparency of the atmospheric channel. The 
power of the transmitted laser beam is ® L . 

From (3) and (4) the expression for p max can be derived: 


P 


max 




TLp^.d-e- 2 )/^ ’ 


(5) 


where e = 2.7182 

The value of the beam radius p z , which allows the extreme 
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magnitude of p max is 

P - = P z,opt 


2.T,.T a .<I>L 


Ti.e./ 


/min can be calculated from the condition 


O 


^min 


P d I SNR=const 


( 6 ) 


(7) 


n.x I .R r 

When SNR value is given, ® pd is calculated from the bellow 
expression [3]: 

( 8 ) 
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C, 


21 B T.A 

Rvu 


+ e“.i?,.(o pd +® B ) 


where R r is the integral sensitivity of the photodetector: 

7? i (^ 0 ) = 8,06.10 5 ti(^ 0 K (9) 

From (8) and physical considerations, the following 
expression for optical signal flux at the entrance of the 
photodetector is reached 
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( 10 ) 


Background optical flux ® B is calculated using the equation 

~ D 2 n2 a a (H) 




Optimal value of the beam divergence angle is 

(U'^fMrad] (12) 

Considering formulas (6) through (12) the final expression 
for 0 tj opt is reached 


0 


t, opt 


1 2T ( T fl T r 7? r O L 

Z V eO PD l SNR=const 


(13) 



-♦-Arithmetic mean -«-Mean + o -A-Mean-a 
Figure 4: Generalized model of the atmospheric visibility 

The equations, which approximate the tree models for P(S M 
> S M , min) illustrated in figure 4, are: 

- worst case “Mean - a”: 

P(S M >S M , min )= 0.00004 x 3 - 


(14) 


; M — u M, min > 

- 0.0019 x 2 - 0.0077 x + 0.9972 

- arithmetic mean: 

P(S M >S UMn ) = 0.0000006c 4 - ( 15 ) 

- 0.0000 lx 3 - 0.0007X 2 - 0.0053X + 0.9991 

- best case “Mean + a”: 

P(S M >S M ,mJ = -0.0005x 2 + 

+ 0.00 16c + 0.9975 
For a complete statistical model, the cumulative 
distribution function ( CDF) and the probability density 
function ( pdf) of S M are needed. By definition CDF is equal 
to the probability P(5 m - *^m, min), which means: 

CDF — P(S M < 5^^ ) = 1 - P(Sm — ^M,min ) 


(16) 


(17) 


Also 


^ S M min ) = j [/ (S Mmin )dS M (18) 


4.2. Statistical Model of the Atmospheric Visibility 

Real yearly measurements of the atmospheric visibility [9] 
were used for creating this generalized version of the 
statistical model of atmospheric visibility. 

The steps for the generalization of the measured data are: 

- Calculating the arithmetic mean value of each P(S M > S M , 

min) 

- Calculate the standard deviation a for the given statistical 
data 

- For each Sm > Sm, min set the interval [x - a; x + a], where 

x = 5 m . 

- Plot the results as a graphic and approximate the resulting 
curves using the simple trend line tool 


The probability density function f(S M ) for each statistical 
model (arithmetic mean, mean + a and mean - a) is equal 
to 


f(S m) 


dP(S M — ^M,min ) 
dS TV/T 


(19) 


Having all the functions P(S M > 5 M ,min), CDF(S U ) and f(S M ) 
defined, the statistical model for the atmospheric visibility 
is completed. In the remainder of this paper the polynomial 
representing the arithmetic mean, (15), will be used as a 
model for P(S M > S M> min). 


4.3. Statistical Model of the Random Vibrations in the 
Initial Direction of the Optical Beam 


23 


In order to define a model for FSO availability, it is also 
needed to have the model of the random jitter in the initial 


Tsvetan Mitsev et al., International Journal of Wireless Communications and Network Technologies, 5(3), April - May 2016, 21-27 


direction of the laser beam axis. 

Consider the case illustrated in Figure 2. Having no 
statistical data in the literature, it is practical to assume a 
Gaussian distribution of the angular misalignments y (y max 
= Qmax) of the optical beam axis from its original direction 
[2, 12]. Taking into account (1), it can be concluded that the 
linear displacements Ap in the plane z = Z are also normally 
distributed: 


/( Ap) = 


1 


V27TG 


-exp 


Ap 


Ap 2 


'Ap J 


Also 


2 _ ryl 2 

a Ap — Z O y 


The cumulative distribution function is 

Pmax 

CDF = F(p lmx )= jf(Ap)dAp-. 
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1 + erf 


Ap 


\ 


a/2c 


Ap 


m 


( 21 ) 


( 22 ) 


4.4. Mathematical Model of FSO System Availability 

In order to have an FSO system working reliably, it is 
required that the power of the optical radiation, at the 
receiver’s site (plane z = Z), is greater than or equal to the 
minimal acceptable power, so 7(p, z) > /min [1]. 

FSO will be in outage if: 

1) . S M falls below a critical value S M , min, according to which 
the system parameters are calculated, 

Or 

2) The linear misalignments Ap exceed the maximal 
acceptable value p max , calculated using (5). 

Condition 1 says that the probability for reliable work of the 
FSO system is P(S M > 5m, min) [1], which was already 
defined in equations (14) through (16). 

Condition 2 means that FSO will work reliably if the 
random linear shifts in the direction of the optical radiation 
are less than p max (Ap < p max ). So the reliable work of FSO, 
depending on Ap is defined by 


max r max 

j f (Ap)dAp = 2 | / (Ap)dAp = 


(23) 


= 2( F( p max ) - F (0)) 

Where E(Ap) is the cumulative distribution function of the 
linear misalignments (22) and p’^is equivalent to I - p max I. 
Considering 1) and 2), Availability is equal to the joint 
probability 

Availabili try = 

(24) 

— ^M,min ’ Pmax ~ ~ Pmax ) 

Knowing that/(S M ) and/(Ap) are independent the above 
equation is equal to [6]: 


P(S M > Vmta/Pmax ^ A P^Pmax) = 

qo Pmax 

= J J f (S M ; Ap) dS M dAp = 


M,min P n 


(25) 


r max 

2 \f{S u )dS u J / (Ap ) dAp 


= 2 P(S M > S M vmin )P( 0 < Ap < p max ) 

The needed system parameters are calculated as defined in 
(4) through (12). 

Using equations (14) through (19) and (23) through (25), a 
simple analytical expression for FSO availability is derived: 
Availabili ty = 

= 2(F(p max )-F(0))P(S M >S MMn ) <26) 

The system reliability calculated using the above formula 
will give the percentage of time, during which the system 
will work correctly using the calculated optimal system 
parameters. 

5. SIMULATION RESULTS 

In order to calculate FSO availability, the following typical 
system parameters were used: channel capacity, C\ - 1 
Gbps; quantum efficiency of the photodetector material 
rj(A, 0 ) = 0,7; SNR = 11,2 (corresponds to BER=10' 8 );central 
wavelength = 1,55 pm; T = 300 K; aperture coefficient A 
= 5; value of the resistor in the feedback of the preamplifier, 
R Fb - 1 kO; x r = x t = 0,85; R r = 5,5 cm; transmission 
wavelength of the interference filter before the 
photodetector A^f = 10 nm; background radiation, L^ B = 
10' 2 W/m2.sr.Ang (corresponds to bright day); angular 
width of the receiving antenna 0 r = 5 mrad 

Figure 5 shows the availability of FSO system using fixed 0 t 
(Figure 5a), and variable beam divergence angle 0 t = 0 t>opt = 
var (Figure 5b). 


0 t = 0,609 mrad = const 



8 <P L , mW 


S\ f, km 


a) 
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b) 

Figure 5: Availability of FSO system using a) fixed beam 
divergence angle and b) optimal beam divergence angle. Sm= var, 
Ol = var, z = 2 km 


The figures 5a) and 5b) represent the percentage of time, in 
which the system will work reliably with optimal system 
parameters calculated according to the variable S M and O l . 
It is observed that there is about 10 to 20 % gain in the 
system availability, if working with optimal value of 0 t . 
Figure 6 a) through Figure 6 c) show a comparison between 
the availability of FSO using 0 t = 0 t , op t = var and a system 
that uses a constant beam divergence angle. Simulations are 
performed for various atmospheric conditions (different 
*5m)- 


It is observed that in all cases using optimal beam 
divergence angle result in higher availability. However the 
results presented in Figure 6 a) through c) also show that 
when using collimated beam constantly (Figure 6 a)) the 
FSO availability is much lower compared to using 0 t = 0 t opt . 
This is because non-divergent optical beams allow minimal 
shifts between optical beam direction and the center of the 
receiving aperture; availability is reduced, because there is 
higher probability of having Ap > pmax. 



a) 



— 0t= 0t,opt 


-0t = const = 0,852 
mrad 


b) 



c) 

Figure 6: Comparison between FSO using different fixed beam 
divergence angles and FSO using 0 t = 0 t , op t = var 


On the other hand, when using larger fixed values of 0 t , 
availability is closer to optimal case (difference is about 
4%), but this is valid only when atmospheric conditions are 
good (Sm > 5 km). When the transparency of the 
atmosphere in the communication channel drops the system 
will not work reliably. 
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®t opt 



b) 

Figure 7: Availability of FSO system using: a) fixed beam 
divergence angle; b) optimal beam divergence angle (Sm= var, z = 
var, O l = 10 mW) 

Figure 7 presents the change in FSO availability under 
different atmospheric conditions and different lengths of 
the communication channel (the distance between 
transmitter and receiver). The results confirm the previous 
findings that using collimated optical beam (with fixed 
beam divergence angle) results in lower FSO availability 
compared to using FSO capable to keep beam divergence 
angle to its optimal value (0 t = 0 t , opt = var). 

Figure 8 shows FSO system availability when different 
receiving antennas are used, calculations of FSO 
availability are performed for 0 t = 0 t , op t and assuming that 
there is a jitter in the direction of the optical beam 
propagation. 



Figure 8: Availability of FSO system using optimal beam 
divergence angle, when different receiving apertures are used (z = 
2 km, O l = 10 mW) 


System parameters are calculated using equations (5) 
through (13). Availability is calculated as defined in (26). 
As expected, when smaller apertures are used system 
availability is lower. This is because using smaller 
receiving aperture requires more optical power at the plane 
of the receiver z = Z, which means that working with 
smaller beam divergence angle 0 t is needed. In this case 
FSO system availability is reduced because there is higher 


probability of having unbearable shifts from the initial 
direction of the optical beam. 

Aperture averaging could be used to further improve the 
system performance (this optimization is not analyzed in 
the paper). Antennas with larger radiuses can also be 
utilized to minimize the impact of the scintillation within 
the optical beam, caused by atmospheric turbulence. 

The outage probability of FSO system can easily be 
calculated as (1 - Availability). This is the amount of time 
the system will work unreliably. 



Figure 9: Outage probability of FSO system using optimal value 
of 0 t and z = 2km, Ol = 10 mW 


Figure 9 depicts the outage probability of FSO system 
working with optimal value of 0 t and in the presence of 
random jitter in the direction of the optical beam 
propagation. Calculations are performed, considering the 
three models of S M ((14) through (16)), presented in Figure 
4. It is observed that, when atmospheric visibility is low and 
system is working with smaller values of 0 t; opt as calculated 
from (13), the outage probability is high, because there is 
higher probability of having Ap > p max . On the other hand, 
when transparency of the atmospheric channel is higher, 
FSO works with wider beam divergence angles. In this case 
there is a higher probability for system outage due to 
unbearable decrease of the atmospheric visibility Sm- 

6. CONCLUSION 

In this paper a generalized statistical model for atmospheric 
visibility was presented. A mathematical model for 
evaluating the FSO systems availability was proposed. Using 
that model the availability of FSOs with different system 
parameters was analyzed. A comparison of the availability of 
FSO system using fixed 0 t and FSO using 0 t = 0 t; op t was 
made. Simulations were performed considering different 
atmospheric conditions. The results clearly show that using 
optimal values of FSO parameters significantly increases the 
availability of the system. 

The presented models can be used to properly choose FSO 
parameters (beam divergence angle, channel length, receiver 
aperture) so that the communication system will work 
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reliably, or they can be used to analyze the availability of a 
commercial FSO system mounted in a given geographic 
region (with known atmospheric conditions). 

A possible future development of this work is the creation of 
an algorithm for real time monitoring of the channel 
parameters, which could be used as a feedback for tuning 
the system parameters to their optimal values by using 
adaptive optical terminals. It is expected that this would 
enable the FSO system to maintain a given BER throughout 
most of the time, even in poor atmospheric conditions. 
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